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Abstract 

In this paper, we consider polynomials orthogonal with respect to a varying 
perturbed Laguerre weight for t < 0 and z on certain contours in 

the complex plane. When the parameters n, t and the degree k are fixed, the Hankel 
determinant for the singular complex weight is shown to be the isomonodromy r- 
function of the Painleve III equation. When the degree k = n, n is large and 
t is close to a critical value, inspired by the study of the Wigner time delay in 
quantum transport, we show that the double scaling asymptotic behaviors of the 
recurrence coefficients and the Hankel determinant are described in terms of a 
Boutroux tronquee solution to the Painleve I equation. Our approach is based on 
the Deift-Zhou nonlinear steepest descent method for Riemann-Hilbert problems. 
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1 Introduction and statement of results 


Let Wt{x) be the following singularly perturbed Laguerre weight 
Wt{x) = w{x] t) = X e (0, +oo) 

with 

t 


(1.1) 


( 1 . 2 ) 


(1.3) 


VAx) = X — log a; H—, t > 0. 

X 

The Hankel determinant is dehned as 
Dk[w{x-,t)] = det(/ii+j)Jjio, 
where fij is the j-th moment of Wt{x), namely, 

roo 

jjij = / x^Wt{x)dx. 

Jo 

Note that when f > 0, the integral in the above formula is convergent so that the Hankel 


determinant Dh[w]t] = Dk[w{x]t)] in (1.3) is well-dehned. Moreover, it is well-known 
that the Hankel determinant can be expressed as 


fc-i 

Dk[w]t] = 

j=0 


(1.4) 


see [26l p.28], where is the leading coefficient of the fc-th order polynomial or¬ 


thonormal with respect to the weight function in (1.1). Or, let iTk^nix) be the /c-th order 
monic orthogonal polynomial, then 'jk,n{t) appears in the following orthogonal relation 

poo 

/ Trk,n{x)x^= 'yk^{t)6jk, j = 0,1, ■ ■ ■ , /c 

Jo 

for hxed n. Moreover, the monic orthogonal polynomials 'Xk,n{x) satisfy a three-term 
recurrence relation as follows: 

X'Xk,n{x) = 7rfc+i,„(x) Q;fc,n(t)vrfc,n(x) -f IJk,n{^)T^k-\,n{x), A) = 0, 1, • • • , (1.5) 

with 7r_i^„(x) = 0 and vro,n()c) = 1, where the appearance of n and t in the coefficients 


indicates their dependence on n and the parameter t in the varying weight (1.1). 


In this paper, however, we will focus on the case when t < 0. Since all the above 
integrals on [0, cx)) become divergent for negative t, we need to deform the integration 
path from the positive real axis to certain curves in the complex plane. Consequently, 
the orthogonality will be converted to the non-Hermitian orthogonality in the complex 
plane. More precisely, let us dehne the following new weight function onT = riUr 2 Ur 3 : 


Wt{z) = w{z]t) = CjO 




, z G r,-, with Cl = 1, C 2 = CK, C 3 = 1 — a. 


( 1 . 6 ) 


2 





where a is a complex constant, the curves Fi = {26, oo), r 2 = {(5(l + e®®)| 9 G (0, vr)} and 
Fa = {h(l + e*®)| 9 G (— 7 r, 0 )}; see Figure 6 being a positive constant. The potential 
is dehned in the cut plane C \ (—oo, 0] as 


Vt{z) = z — logz +-, aigz e {—7i,7i), t<0. 
z 

The orthogonality relation now takes the form 


Trk,n(:):^w,{z)dz = j =0,1,--- ,k. 


( 1 . 7 ) 


( 1 . 8 ) 



Figure 1: Contour of orthogonality, F = FiUF 2 UF 3 . 


With the weight function Wt{z) given in ( 1.6[ ), the corresponding Hankel determinant 
Dn[w; t] in (1.3) is well-dehned. However, since Wt{z) is not positive on F, the orthogonal 


polynomials TTk,n{z) in (1.8) may not exist for some k, and (1.4) only makes sense if all 
polynomials iij^n for j = 0,1, • • • ,k — 1 exist. It is worth mentioning that as part of our 
results, we will show that there exists a < 0 , such that 7in,n{z) exists for n large enough 


and t > ter', cf. Section 2.1 The recurrence relation (1.5) still makes sense for such t if 


all of Tik_i^n{x), T^k,n{x) and 'Kk+i,n{x) exist. Note that in the literature, the polynomials 
with non-Hermitian orthogonality have been studied in several different contexts; see for 
example HEIEIEI, where the cubic and quartic potentials are considered. 

One of the main motivations of this paper comes from the Wigner time-delay in the 
study of quantum mechanical scattering problem. To describe the electronic transport 
in mesoscopic (coherent) conductors, Wigner [22] introduced the so-called time-delay 
matrix Q] see also Eisenbud [15] and Smith [25|. The eigenvalues Tk of Q, called the 
proper delay times, are used to describe the time-dependence of a scattering process. 
The joint distribution of the inverse proper delay time 7 ^ = I/ta, was found, by Brouwer 
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et al. [8], to be 


^( 71 , 72 , •••7n) = n (1-9) 

^ i=l 

Then the probability density function of the average of the proper time delay, namely 
the Wigner time-delay distribution, is defined as 

-I«n / n ^ \ n 

a(r) = ^/ n n‘*T‘- (!■“) 

n Jwi l<i<j<n \ 1 9'*y 


The moment generating function is the Laplace transformation of the Wigner time-delay 
distribution 


Mn{z) = 


nOQ 1 /* 

/ e-"PjT)dT^— '[{ite-"-'--" n h. 

Jn i=l l<i<j<n 


which is closely related to the Hankel determinant (1.3) as follows: 


( 1 . 11 ) 

i=l 


Mn{nt) 


Dn[w{x]t)] 


( 1 . 12 ) 


Recently, Texier and Majumdar IZ71 studied the Wigner time-delay distribution by using 
a Coulomb gas method. They showed that 


^^(7i 172 , ■ ■ ■ , 7ii) exp{ -v?E\p{x )]} for large n, 


(1.13) 


where p{x)dx is the unique minimizer for an energy problem with the external held Vt{x) 
in (1.2), and E[p{x)] is the minimum energy. Moreover, the density p{x) is computed 
explicitly in [27], namely. 


p{x) = ^ y/{x — a){b — x) for x E [a,b], with 0 < a < 6, c = t/V^. (1-14) 


2tix‘^ 


Here positive a and b are independent of x and implicitly determined by t as follows: 

= 3. (1.15) 




One may notice that p{x)dx is a probability measure on [a, b] as long as a -|- c is non¬ 
negative. Since a -|- c is a continuous function of f, we see that p{x) in (1.14) is non¬ 
negative for t > ter-, where ter = — \ (2^/^ — l) is the critical value of t corresponding 
to the case a -f c = 0; see Theorem]^ It is very interesting to observe that, for this 
fer < 0, we have Ccr = — Ocr < 0 and 


p{x) = 7 ^—- acrf^‘^{bcr - 


2tix‘^ 


(1.16) 
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where a phase transition emerges at the left endpoint x = Qc 
and her are explicitly given in (1.32) and (1.33). 


Here the critical valnes 


t. 


a. 


cr 


It is also interesting to look at onr problem from another point of view. Dne to the 

X 


term - in the exponent of (1.1), we may consider the origin as an essential singnlar 


point of the weight fnnetion. In recent years, orthogonal polynomials whose weights 
possess essential singnlarities have been stndied extensively. For example, Chen and Its 
[9] consider orthogonal polynomials associated with the weight 


Wt{x) = x°‘e 


X G (0, cxd), a > 0 and t > 0. 


(1.17) 


They show that, for hxed degree n, the recnrrence coefficient satishes a particular 
Painleve III equation with respect to the parameter t, and the Hankel determinant of 
hxed size Dn[wt{x)] equals to the isomonodromy r-function of the Painleve III equation 
with parameters depending on n. The matrix model and Hankel determinants Dn[w,t] 


associated with the weight in (1.17) were also encountered by Osipov and Kanzieper 


in bosonic replica held theories. Later, the large n asymptotics of the Hankel determi¬ 


nants Dn[wt{x)] associated with the weight function in (1.17) is studied by the current 
authors in isi and For t G (0,d], the asymptotics of the Hankel determinants are 
derived and expressed in terms of certain Painleve HI transcendents. The asymptotics of 
the recurrence coefficients are also obtained therein. In the case of the Gaussian weight 
perturbed by essential singularity 




X G 


(1,18) 


the double scaling limit of the Hankel determinants are also characterized in terms of 
Painleve HI transcendents by Brightmore et ah in [7]. Recently, Atkin, Claeys and 
Mezzadri [1] extend the results to the case of Laguerre and Gaussian weight perturbed 
by a pole of higher order at the origin, they obtain the double scaling asymptotics of the 
Hankel determinants in terms of a hierarchy of higher order analogs to the Painleve HI 
equation. 

The main objective of this paper is to study the Hankel determinant Dk[w]t] with 


respect to the weight (1.6) in the region t < 0. First, for hxed degree fc, we will show that 
the recurrence coefficient ak^n satishes a Painleve HI equation, and the Hankel determi¬ 
nant Dk[w] t] equals to the isomonodromy r-function of the Painleve HI equation. Then, 
we will derive the double scaling limit of the Hankel determinant Dn[w]t], the recurrence 
coefficients and leading coefficients of the associated orthogonal polynomials. Our results 
are described in terms of a certain tronquee solution of the Painleve I equation. 


1.1 A model Riemann-Hilbert problem for Painleve I 

To state our results, we need certain special solutions to the Painleve I equation 

2/"(s) = 6|/2(s)-1-s. (1.19) 
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The reader is referred to m Ch.32] for properties of the Painleve I equation, as well as 
the other Painleve equations. In [21] , Kapaev formulates the following model Riemann- 
Hilbert (RH, for short) problem for T((^) = associated with the Painleve I 

equation. This model RH problem will play a crucial role later in the construction of a 
local parametrix in the steepest descent analysis. 

(a) T((^; s) is analytic for ^ G C \ P^, where 

r^r = 7-2 u 7_1 U 7i U 72 U 7* (1-20) 

are illustrated in Figure 



Figure 2: The contour Fij, associated with the Painleve I equation 


(b) Let denote the limiting values of as ( tends to the contour F^ 

from the left and right sides, respectively. Then, \h((';s) satishes the following 
jump conditions 


T+(C;s) = d/_(C;s)<^ 


1 Si 
0 1 

1 0 
i 1 
0 -i 
-i 0 


2 ; e 7fc, k = ±1; 
2 ; e 7fc, k = ±2; 
2 ; e 7*, 


where si = (1 — a)i and s_i = m, with a being a complex constant, 

(c) As C —^ cxD, tI/(C;s) satishes the asymptotic condition 

3^'{C:S)cr3 




72 


7C 


c 


3^ 

21 I e“ 


( 1 , 21 ) 


( 1 . 22 ) 
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for argC G (—7r,7r), where 


0{C,s) = +sCK 
(Ti and (73 are the Pauli matrices 


ai = 




(1.23) 


It is known that, for each a E C, 
Vais) = 2('h_2(s))i2 


(1.24) 


is a solution of the Painleve I equation (1.19). As a consequence, the above RH problem 
for s) has a solution if and only if s is not a pole of ya{s). Due to the meromorphic 
property of the Painleve I transcendents, one also see that the solution of the above RH 
problem for s) is meromorphic in the parameter s. Moreover, it is shown in Kapaev 
IZH that ya{s) is the so-called tronquee solution of Painleve I whose asymptotic behavior 
is given by 


ai 


ya{z) = yo{z) -^2 8 (- 32 ;) 8 gxp 


TT 


1 11 1 , ,5 

--2 4 34 (- 2;)4 

5 


(^l + 0(z-i)) (1.25) 


as z —>■ oo and aigz = arg(—z) +-n E [|7r,7r]. Here yo{z) is the tritronquee solution 


satisfying 

yo{z) ~ \J-zl<o 


1 -|- afc(—^ 


\-bkl2 


k=l 


TT TtT 

as 2 ; —)■ cxD,-< arg 2 ; < —, (1-26) 

5 5 


where the coefficients can be determined recursively; see for example Joshi and Kitaev 
1201. The solution ya^z) will appear in our main results below. 

We mention several known facts about the coefficients in (1.22) in addition to (1.24). 
For example, the explicit formulas of \h_i(s) and \['_ 2 (s) are given in [21] as 


T_i(s) = -'R„(s)cr 3 , 


'I'-2(s) = § I + yo:(s)ai |, 


(1.27) 


where 

= \y'a{s) - 2yl{s) - sya{s) (1.28) 

is the Hamiltonian of Painleve I. 
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1.2 Statement of main results 


First of all, when the degree k is hxed, we show that the recurrence coefficient ak^n{t) 
satishes a particular Painleve III equation with certain initial conditions. Moreover, we 
prove that the Hankel determinant Dk[w{z] t)] is related to the r-function of the Painleve 
III equation. Similar results for the weight in (1.17) have been obtained by Chen and 
Its [9]. 


Theorem 1. For fixed non-negative integer k, let ak^n the recurrence coefficient in 

2/c + 1 + n 


(1.5), and 

^k,n{t) (^k,n{t) 


n 


Then ak{t) = Okffit) satisfies the following Painleve III eguation 


// 

— 


AY A /nt 1 

+ n{2k + ^ + n)^ + ^ + — -, 

Ok t F F t ttk 


with the initial conditions 0^(0) = 0, a'^(O) = 1. Moreover, we have 

(] = const ■ t(«) 


(1.29) 


(1.30) 


(1.31) 


where T{t) is the Jimbo-Miwa-Ueno isomonodromy r-function of the above Painleve III 
eguation. 

Next, we let A: = n and consider the double scaling limit when n ^ oo and t —?• ter 
simultaneously. We show that the asymptotics of the Hankel determinant Dn[w{z;t)] 
associated with the weight in (1.8) can be expressed in terms of the tronquee solution 
yffis) of Painleve I equation given in (1.24). 

Theorem 2. Let the constants ter, cter o.nd ber be defined as 


= --(2^/3 _ _ _o.051 


and 


a„ = 1(3 - 2'-'® - 2=''“) cs 0.076, b„ = 1(1 + 2'-'“ + s. 5 771 ^ 


(1.32) 


(1.33) 


For n —>■ 00 and t —)■ ter in a way such that 

S (flOer^ ^ {perber^ ^ ificr ^cr') ^ ^ U ^ (fer A 

remains bounded. Suppose a & C is fixed and s* is not a pole of the tronguee solution 
yffis), then an asymptotic approximation of the logarithmic derivative of the Hankel 
determinant Hn,n = Dn[w;f\ associated with the weight function (1.6) is given by 


d ( Her I ber ^ 1 2(6c 


ns (2aer) ^'/affib'e 


2/„(s*) + 0 — . (1.35) 


ns 












We would also derive the double scaling limit of the recurrence coefficients and the 
leading coefficients of the orthonormal polynomials. 

Theorem 3. Under the same conditions as in the previous theorem, the monic polyno¬ 


mial TTn,n{z) defined in (1.8) exists for large enough n and t close to ter- Moreover, we 


have the asymptotics of the recurrence coefficients 


and 




fin,n 


^2 = 
/ n,n 


_ 1 

\ Qjcr^cr \ 

{per 

(Zer) 


16 


2 


W 


, /c — \~ o 

— OcrYUn^/^ 

{flacrificr - acr)Y^^ya{s*) 1 
4 n2/5 


3/5 


+ 0 


„ 3/5 ) 


—nl 


*^cr) 


1 + 


2{2a,rY/^'H^{s* 


{be 


„)h5nV5 


+ 0 


n 


2/5 


(1.36) 


(1.37) 


(1.38) 


where 'Ha{s) is the Hamiltonian of Painleve I given in (1.28) 


Remark 1. It is well-known that the tronquee solutions of Painleve I are meromorphic 
functions and possess inhnitely many poles in the complex plane. Therefore, to make the 


results valid in the above theorems, we require the s* in (1.34) is bounded away from the 
poles of ya{s). Recently, through a more delicate triple scaling limit, Bertola and Tovbis 
[1] successfully obtain the asymptotics near the poles of ya{s). Similar results near the 
poles of ya{s) might be derived by using their ideas in jl]. However, we do not pursue 
that part. Instead, we focus on the main task of the present paper to demonstrate that 


the Painleve I asymptotics can also occur for the weight (1.1) with negative t 


The rest of the paper is arranged as follows. In Section we provide a RH problem 
for the orthogonal polynomials with respect to the weight ( |1.6 ). A transformed version 
of the solution is shown to fulhll a Lax pair, which is closely related to the Painleve 
III equation. Several differential identities are stated and justihed. Theorem [T] is also 
proved in this section. Section is devoted to the determination of equilibrium measures, 
involving a positive measure and a signed measure. In Section]^ we carry out a nonlinear 
steepest descent analysis of the RH problem for the orthogonal polynomials. Particular 
attention will be paid to the construction of the local parametrix at the critical endpoint 
z = Ocr, where the Painleve I transcendents are involved. Then, the proofs of Theorems 
and are given in the last section. Section 


2 Finite Hankel determinants and Painleve III equa¬ 
tion 

In this section, we state the RH problem for the perturbed Laguerre orthogonal poly¬ 
nomials. Then we show that after some elementary transformations, the RH problem 
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is transformed into a RH problem for the Painleve III eqnation. As a conseqnence, we 
derive a Painleve III eqnation satisfied by the recnrrence coefficient ak,n up to a trans¬ 
lation, and establish a relation between the hnite Hankel determinant of the pertnrbed 


Lagnerre weight in (1.8) with the r-fnnction of this Painleve III eqnation. Several dif¬ 


ferential identities for the Hankel determinants and the recnrrence coefficients of the 
pertnrbed Lagnerre orthogonal polynomials are also derived. The identities are impor¬ 
tant in the asymptotic analysis in later sections. Althongh onr calcnlations are similar 
to those in Chen and Its [5] , we think it is convenient for the reader to have more details. 

2.1 Riemann-Hilbert problem for orthogonal polynomials and 
differential identities 

We state the RH problem for the pertnrbed Lagnerre orthogonal polynomials as follows: 
(Yl) Y(z) is analytic in C\Pj, j = 1,2,3; see Fignre[^ 

(Y2) Y{z) satishes the jnmp condition 

1 w{z;t) 


Y+{z) = YZz) 


0 


j e r = Ti U Fj U Fa 


( 2 . 1 ) 


where w{z]t) is the weight fnnction piecewise-dehned in (1.6); 
(Y3) The asymptotic behavior of Y{z) at infinity is 

z^ 0 


Yiz) = iI + Oil/z)) 


0 z 


-k 


as z ^ oo] 


( 2 . 2 ) 


(Y4) As j ^ 0, Y(z) = 0(1). 

Using a by now standard argnment, originally dne to Fokas, Its, and Kitaev im, the 
solntion of the above RH problem, if it exists, is nniqnely given by 


Y(z) = 


T^k[Z) 


1 r 

27Ti Jr S-2 


-2iri72_i irt_i(z) - 7^.1 / 


(2.3) 


where 7ik{z) = 7rk,n{z) is the monic pertnrbed Lagnerre orthogonal polynomials dehned 


in ( 1.8[ ) and 'jk = 'Ik,nit) is the leading coefficient for the ortho normal polynomial of 
degree k. To show the existence of 7Tn,niz) when n is large enongh, we will apply a 
series of invertible transformations to transform the original RH problem Y to a new RH 


problem for R, which is solvable for snfhciently large n, t close to t^r as in (1.34), and s 


is not a pole of the tronqnee solntion |/a(s). Tracing back the invertible transformations, 
we will see that the RH problem is solvable nnder the same conditions. Indeed, it is also 
possible to prove the solvability for f > fcr- However, since we are interested in the phase 
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transition near ter, we don’t consider the case when ter < t < 0 in the snbseqnent analysis. 
Thns, the pertnrbed Lagnerre orthogonal polynomials are well-dehned for k = n large 
enongh. 

Next, we derive some differential identities for the recnrrence coefficients and the 
logarithmic derivative of the Hankel determinant associated with the pertnrbed Lagnerre 
weight w{z) = w{z;t) in (1.6). The resnlts are expressed in terms of the entries of Y{z). 


Lemma 1. Assume that f > 0. Let ak,n{t) and /3k,n{t) be the recurrence coefficients in 

(2.4) 


(1.5), and Dk[w,f\ he the Hankel determinant in (1.3). Define 


and 


= afc,n(t)-( 2 fc + 1 + n) 

n 


Hk,nit) = t^^\ogDk[w]f\. 


Then we have 


ak,nif) ^'hk,r 




dz = 2witjlJt)Yn(0)Y\2(0), 


(2.5) 


( 2 . 6 ) 






k{k + n)+ t^Hk,nit) - Hk,nit) 


(2.7) 


and 


—Hk,n{t) = n / --- dz = -n yi 2 ( 0 )F 2 i( 0 ). (2.8) 


d 

dt 


Proof. Since f > 0, the orthogonal polynomials TTk^n exist for all nonnegative k and 
positive n. First, we consider the recnrrence coefficient ak,n{t). Based on the three-term 


recnrrence relation (1.5) and the orthogonality condition (1.8), we get 


ak,n{t) = 'ylnit) / Z7rl^{z)w{z)dz. 


(2.9) 


Using the fact that w{z) = -f and integrating by part once, the above 

formnla gives ns 


akAt) = tilAt) / dz. 

Jv ^ 


( 2 . 10 ) 


Then (2.6) follows from a partial fraction decomposition of orthogonality 


condition (1.8), and the explicit expression of Y{z) in (2.3). 
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Next, we consider the Hankel determinant. Recall that the Hankel determinant can 
be expressed in terms of the leading coefficients as 

fc-i 

Dk[w]t] = 

j=0 


see (1.4). Taking logarithmic derivative of both sides of the above equation with respect 
to t and using the integral representation of the leading coefficients in (1.8), we get 

fc-i 

Hk,n{t) = (2.11) 

j=0 


Differentiating the above formula again, we get from (2.4) 

fc-i 

dt ' dt 


j=0 


Let pk,n(t) be the coefficient of the ^ term in 7rk,n{z), i.e., 
TTkA^) = + Pk,ni^)z’'~^ H-. 


Comparing the powers in the recurrence relation (1.5), we obtain 
(^k,n Pfc,n(^) Pfc+l,n(^)- 


( 2 . 12 ) 


(2.13) 


(2.14) 


To derive one can see from (2.12) and (2.14) that it is sufficient to obtain 

^Pk,n{t)- This can be done by taking derivative of the following orthogonal formula with 
respect to the parameter t 

/ 7^k,n{^)'^k-l,n{^)'^i^)d^ = 0 - 


More precisely, taking into account the orthogonal relation (1.8) and the fact that 
= -riyj{z;t), we have 


^pfc.n(t) = 


T^k,n{^)T^k-l,n{^)w{z] 


dz. 


(2.15) 


Then, (2.8) follows from a combination of (2.12), (2.14) and (2.15), as well as the deh- 
nition ofF(^) in (2.3). 


Finally, let us study Using the ideas leading to (2.10), we have 

'7lk,n{^)T^k-l,n{^)'w{z) 


Pk,nit) t'Jk-ir 


dz -Pfc,n(^), 


n 


(2.16) 


where pk,n{t) is introduced in (2.13). The hrst term on the right-hand side is 

cf. (2.8). An expression for the term on the extreme right can be obtained by deriving 

Pk,n{^) = from (|2.14[), and using ([^ and (|2.11[). 


This completes the proof of our lemma. 


□ 
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Remark 2 . For later use, we need the differential identities of Lemma [T] in the case when 
k = n is large and t ~ ter- They can be obtained through an analytic continuation 
argument. Indeed, Y{z) determined by RH problem exists in this case, and is related 
to the T-function of the third Painleve equation after an elementary transformation 


given in (2.17). Thus Y{z) is meromorphic with respect to t in the cut plane arg(—t) G 


(— 7 r/ 2 , 37 r/ 2 ). In particular, both Y and the Hankel determinant are analytic in a domain 
containing the interval t > 0 and a neighborhood of t = ter- Note that the identities 


( 2 . 6 )-( 2 . 8 ) hold for t > 0 , then, by analytic continuation, they also hold for t close to 


ter- We conclude that for k = n large and t ~ ter, the identities (2.6)-(2.8) are also true. 


Similar argument has previously been used in Bleher and Deano [5]. 


2.2 Relation to the Painleve III equation 

Introduce a purely imaginary parameter s = niy/Yt^ and define 

/o\\ 


F — — , A ^r, 

ni J \ni 


(2.17) 


where T* = T^ U T^ U Tg = :^r is the rescaled contour. Then, <I>(A) = d’(A,s) solves 
the following RH problem with constant jumps: 

(i) ‘h(A) is analytic for A G C \ U^^gT*. As T* and T only differ by a scale, one may 
refer to Figure to see the properties of the contour F*. 


(ii) ‘h(A) satisfies the jump condition 


$+(A) = $_(A) 


1 9 

0 1 


, A G F* j = 1,2,3, 


(2.18) 


where Ci = 1 , C 2 = a, C 3 = 1 — a; cf. ( 1 . 6 ) 


(iii) The asymptotic behavior of *h(A) at inhnity is 


4>(A)= L + X] 


k=l 




A(t+^) 


t73g —0-3 


as A —)■ cx). 


where 


27ri7fe,„(p 


(2.19) 



-(t+fcVs 


, ( 2 . 20 ) 


In the above formula, and pk,n are, respectively, the leading coefficient of the 
fc-th orthonormal polynomial, and the coefficient of the z^~^ term in the /c-th monic 


orthogonal polynomial introduced in (2.13), with respect to the varying perturbed 


Laguerre weight in (1.6) and (1.8). 
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(iv) The asymptotic behavior of *h(A) at A = 0 is 

$(A) = $(0) + as A ^ 0, 

where 


( 2 . 21 ) 


k=l 


m = ri) 


j\ (f+^Vs 


Qk,n{t') 1 C}i,n{^')(ik,n{i^ ' \ S 


-^nas 


, ( 2 . 22 ) 


with 


/,x 7rfc(0) f 7rk{z)w{z]t) o ^ /,^7rfc-i(0) 

Ck,n{t) = = ^T^^lk-l,n\t) -77^- 

Itti Jy z 7rfc((J) 

Now, from the above RH problem, we derive the following Lax pair for the function 
$(A, s), which is exactly the same as the Lax pair for Painleve III; see [TSl pp. 195-203]. 


Proposition 1. For the matrix function <h(A,s) given in (2.17), we have 
$a = R(A,s)<1) and ^s = B{X,s)^, 

where 

^-1 ^-2 , TA B_i 

R(A, s) — ~^<^3 H— ^ —1“ o-nd B{X, s) — ~^<^3 + Bq -\— 

Here, the coefficients in the above formula are given below 

. n-\-2k 

/ ^ + K I — 

71-1 = 


f+ fc 




(f)’ 


^nm^Li,n (f) 


—n—2k 


.H — Jc 
2 ^ 


(2.23) 


(2.24) 


(2.25) 


n-\-2k 


and 


A I 1 nQk,n 2Cfc ( ) 

R-2 = — / -x -n-n 

2 \ -2gfc,n(l - Ck,nqk,n) (f) 2Ck,nqk,n “ 1 


Bo - 1 (a-i - (I +1;) ffs). 


(2.26) 


(2.27) 


Proof. Note that the jump matrices in (2.18) are independent of A and s. This implies 
that both 

yl(A, s) = and B{X,s) = ^s^~^ (2.28) 

are analytic functions of A with only possible isolated singularities at the origin and at 


inhnity. Using the asymptotic expansions in (2.19)-(2.22), we hnd that 


A-i — /c^ (T 3 -|- —[<l)_i, CTs], 24_2 — —<h( 0 )cr 3 $( 0 ) ^ 


(2.29) 
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and 


-Bo — 2 [*^-1) *^ 3 ]) B_i — -—, (2.30) 

where [^, = XY — FX is the commutator. Then direct computations give us the 

results. □ 

It is known in several circumstances that the Hankel determinants admit an interpre¬ 
tation as the Jimbo-Miwa-Ueno isomonodromic r-function for the rank 2 linear system of 
differential equations; see [I6] for the Hankel determinants associated with the exponen¬ 
tial weight and [2113] for more general semi-classical weights. Now we have established 
the relation between the perturbed Laguerre orthogonal polynomials and the Lax pair 
for the Painleve III equation. Naturally, the associated Hankel determinant is also ex¬ 
pected to relate to the r-function of the Painleve HI equation. Thus we are in a position 
to prove our first result for hxed degree k. 

Proof of Theorem^ According to Proposition]^ $(A,s) satishes the same Lax pair as 
Painleve HI. Then, applying an argument in [T8l (5.3.4),(5.3.7)], we see that the function 


u{s) — —i{A- 1 ) 12 /{A- 2 ) 12 — 


n 


2mSCk,nll^n 

solves the Painleve HI equation 
{u'f u' 4 


u"{s) = 


u 


-h -(©oM^ + 1 - ©oo) + 


u 


with the parameters ©0 = n and ©00 = —{2k + n). By (2.6), we have 
nt 


u{s) = 


S^k,r 




(2.31) 


(2.32) 


(2.33) 


Substituting (2.33) into (2.32) gives us (1.30). 

Next, we consider the Hankel determinant Dk[w]t]. Denote by ‘hoo('^) and ‘ho(-^) the 
series in the expansions ( |2.19 ) and (2.21), namely. 


$oo(A) = / + 5 ^^ and $o(A) = / + 5 ^<hfcA^ 


k=l 


k=l 


with <l)_i qiven in ( 2 . 20 ) and 


4*1 = f -n^(^k,n{t) - {k^ + nk) + y (1 - 2ck,n{t)qk,n{t)) ) 0-3 + 


IS 


0 * 
* 0 


; (2.34) 


cf. (2.3) and (2.17), where * denotes the off-diagonal entries independent of A. By the 
general theory of Jimbo-Miwa-Ueno [T9|, the isomonodromy r-function for the Lax pair 
in (2.23)-(2.27) is dehned by the formula 

^.-l.^^^^o(A) 


dlogr(s) :=-ResTr<i<ho-i(A)^^dTo(A) -ResTr $-1(A)^1 A|^c?Too(A) [>, 


A=0 [ 


d\ 


A=oo 


a<l)oo(A) 


d\ 


15 






















(2.35) 


where 


% %\ 
dTo{\) = -—asds, dT^{X) = —a^ds; 


see m Eq.(1.23)]. Substituting the dehnition of d>oo(A) and $o(A) into (|2.35D, we obtain 


^ ^Tr($icr 3 -$_icr 3 ) = ^ {-2in?I3k,n + ‘2s‘^Ck,nqk,n - + nk)) . (2.36) 


Now a combination of (2.8), (2.11), (2.14) and (2.16) gives 
^ ldk,n ^ tCk^nQk,n Hk,n 4“ k(^k u). 


(2.37) 


see (2.26) for the definition of Ck,n{t) and qk,n{t)- Thus, we obtain from (2.36) and (2.37) 
that 




(2.38) 


Here use has been made of the relation = n‘^t. In view of the formula (2.5), and 


integrating both sides of (2.38), we arrive at the following relation between the Hankel 


determinant Dk[w,t] and the r-function of the Painleve III equation: 
Dk[w;t] = const • 

which is (1.31[). This completes the proof of Theorem]^ 


□ 


3 Equilibrium measures 


The equilibrium measure with the external held Vt{x) in (1.2) is given recently in Texier 
and Majumdar 


To obtain a double scaling limit at the critical time, we need a mod- 
ihed equilibrium problem, which will involve a signed measure. This signed measure will 
be used to construct the important gf-function and (/)-function in the Riemann-Hilbert 
analysis. The idea of considering a modihed equilibrium problem has been successfully 
applied to study similar double scaling limits in several different problems, such as vary¬ 
ing quartic potentials by Claeys and Kuijlaars [10] and Duits and Kuijlaars [Tl|, and a 
cubic potential by Bleher and Deano [Bj. 


In this section, we will go back to the weight (1.1), consider a regular equilibrium prob¬ 


lem hrst, and see how the critical time occurs. Then, to facilitate our future Riemann- 
Hilbert analysis near the critical time, we will consider a modihed equilibrium problem 
by hxing the left endpoint. This will give us the signed measure we need. 
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3.1 Equilibrium measure and a critical case 

Consider the extremal problem minimizing the energy with the external field Vt{x) in 

(pl): 


r*cxD roo 


I{v) = / Vt{x)di'{x)-\- 


log 


x-y\ 


dh'{x)dh'{y). 


(3.1) 


According to the general potential theory 123!. there exists a unique minimizer dut of 
I{u) among all Borel probability measures du on [0,+cx)), such a probability measure 
is called the equilibrium measure. For the potential Vt{x) = x — log a; + t/x in (1.2), the 
equilibrium measure du^ can be computed explicitly. 

The equilibrium measures for t > ter and t —)■ ter have been computed explicitly in 
Texier and Majumdar |27|. To make the present paper self-contained, we sketch the 
proof below, which differs from that in |2^. Inspired by [2^, and in view of the measure 
for the positive-f case, we assume that the support of dut^x) has only one piece. Also, for 
fixed t, the behavior of the density Vt{x) is expected to demonstrate a weak singularity 
at the endpoints since the contour is deformed to keep away from the possible singularity 
at the origin. We derive the equilibrium measure by solving a scalar RH problem, based 
on the Euler-Lagrange equation (3.4). 

Proposition 2. Let Vt{x) be the density function of the equilibrium measure dut sup¬ 
ported on an interval {a,b) C [0, cxd), such that dutix) = Vt{x)dx. Then fort > ter we 
have 

X c 


Vt{x) = ~ ~ xe{a,b), 


(3.2) 


where c = t/y/ab and a, b are determined by (1.15). Moreover, when t = ter as in (1.32), 
we have 


"(^) ^ ^cr')^ iper x '), X G (^aeriberp 


(3.3) 


where the critical endpoints Ocr and ber are given in (1.33). 


Lagrange equation 


Proof. From (|3.1|), it is known that the equilibrium measure dut satisfies the Euler- 

(3.4) 


Vt{x) + 2 log 


|x - y\ 


diyt{y) = l, xe(a,6). 


where I is the Lagrange multiplier. Differentiating with respect to x, we get 


l//(x)-2p.v. 


My) 


dy = 0, X E (a, b). 


(3.5) 


x-y 

where the integral is taken as the Cauchy principle value. This is an integral equation 
for the density function Vt{x), which can be solved explicitly. Indeed, one can define 

vt{y) 


G{z) = - f 

TT* Ja 


y-z 


-dy, z G C \ [a, 6]. 


(3.6) 
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Then it follows from the Plemelj formula that 


1 /■' 

G'±(x) = —p.v. / 
TT* Ja 


My) 

y-x 


dy ± Vt{x), X G (a, 6), 


(3.7) 


where the integral is the Cauchy principle value. It is readily verihed that G{z) satishes 
the following scalar Riemann-Hilbert problem: 

(i) G{z) is analytic for z G C \ [a, h] , having at most weak singularities ai z = a, 6; 

(ii) G+(x) + G_(x) = -^Vl{x) for x G (a, 6); 

(iii) G{z) ~ as 2 ; —)■ cxD. 

Solving this RH problem yields 

V^M) cz + c 


G{z) = 


27ri 




(3.8) 


where V^{z) = ^ \ c = and c = [^ + 2 (a i)] • attention should 

be paid to the fact that G{z) is analytic outside the interval [a,b], especially at 2 ; = 0. 


Now expanding (3.8) in powers of 1/z, the large- 2 ; behavior of G{z) ensures that 


c = 1 and c-= —3, 


which are indeed (1.15). Furthermore, a combination of (3.7) and (3.8) yields (3.2). 


Moreover, in the critical case when 


Cr.r — 


a straightforward computation gives us (1.32)-(1.33). 


□ 


Remark 3. Of course, the formulas for the density function and endpoints in (3.2) and 


(1.15) hold when t > 0. For any t > 0, the density function is supported on [a, 6] 
with 0 < a < 6 and vanishes like square roots at both endpoints. As a consequence, 
one will obtain usual Airy-type and sine-type asymptotic expansions for the orthogonal 
polynomials near the endpoints and inside the support, respectively. 


3.2 Signed equilibrium measure 

The critical case when t = ter is termed a freezing transition] see 123. One can see 


that, when t = ter, the density function Vt{x) in (3.3) vanishes like a 3/2 root at a. 


This suggests that the local behavior of the orthogonal polynomials near is described 
in terms of the Painleve I transcendents; see [SI Hi]- To precisely construct a local 
parametrix near the endpoint Oer by using the Painleve I transcendents, a delicate study 
near Oer is needed in our subsequent nonlinear steepest descent analysis for the RH 
problem. Therefore, technically it is more convenient to have a measure whose left 
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endpoint of the support is exactly located at acr- Note that in the case when only positive 


measures are involved as in Section 3.1 both endpoints a and b in (1.15) vary when 


the value of parameter t changes. So we need to minimize the same energy functional 


(3.1) among signed measures on [ocr,+ cxd) which are nonnegative except possibly on 
[ocr, dcr + (5i] for some sufficiently small (5i > 0. As there is no symmetry as in la. the 
right endpoint b may depend on t. A similar treatment is also employed in [6]. 

By a similar argument performed in Section |3TT| we hnd the new minimizer explicitly. 

Proposition 3. Let 'iptix) be the signed density function of the minimizer of the mini¬ 


mizer of the energy functional in (3.1) on [acr,+oo). Then we have 
1 




2'kx‘^ 


X 


X — Or 


{x^ + dix + do) , X E [ucr, b], 


where 



do = —t\ —— and di = — 1 



2ar 


t 

^2b 


and b is determined by the eguation 
t 



OtCr j 


2ac_ 


t \ b — c 


= 3. 


(3.9) 


(3.10) 


(3.11) 


It is worth noting that for t = ter, the density of the modihed equilibrium measure 


fjt^x) is reduced to Vcr{x) in ( |3.3[ ). Moreover, near the critical time t = ter, we have 

(brr - aer)~^it - ter) + O {{t - ter)^) , 


b — “t“ 


do = ai - 


-^cr 




(3.12) 


dl = -2aer + \ \ ^{^cr “ ^cr) “ Ar) + O {{t - ter)^) ■ 


Based on the signed measure obtained above, we dehne several auxiliary functions 
which will be used in our further analysis. 


Definition 1. The gf-function is dehned as 

rb 

g{z) = / log(z — s)'0t(s)ds for G C \ (—oo, &], 


(3.13) 


where arg( 2 : — s) G (— tt, tt), and the equilibrium density function ipt{x) is given in (3.9). 
Definition 2. We also dehne the following 0-functions 


2 


^cr[^) 2 


-|- dis do){s — 6)^/^ 
s2(s - 

(S - acr.)^/^(s - 


ds, z E C \ {( — CXD, 0) U (Ocr, cxd)}, (3.14) 


ds, G C \ {(—oo, 0) U (Ocr, cxd)}, (3.15) 


where the branches are chosen such that arg(s — Oer) G (0, 27r), arg(s — ber) G (0, 27r) and 
arg(s — 6) G (0, 27r). 
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From the above definitions, it is immediately seen that g{z) satisfies the Euler- 
Lagrange equation 


g+{x) + g_{x) - Vt(x) -1 = 0, x G {a^r, b), 
and the variational inequality 

‘igiyx) — Vt{x) — I < 0, X E {b, oo), 


(3.16) 


(3.17) 


where I is the Lagrange multiplier introduced in (3.4). Moreover, the ^f-function and 
the (/)-function are related by 


g+{x) - g_{x) = 27ii-2 (0*) (x), x e (a^r, b). 


(3.18) 


Note that (j)t{z) and (j)cr{z) are close to each other when t approaches ter- If we rewrite 
as 


0 cr('^) (t fcr) 0 o('^)) 


then, in view of (3.14)-(3.15), we have 


0o(^) = - acr(l + ro{z)), aTg{z - a^r) E (0, 27r), 

ZCtcf'y dcr^cr 


(3.19) 


(3.20) 


where tq^z) is analytic in a neighborhood of 2 ; = Ocr and ro(acr) = 0. We also need some 
local information of the functions (j)cr{z) and (f>t{z) at critical points 2 ; = Qcr and 2 ; = 0 . 


From their definitions in (3.14) and (3.15), we have 


, / \ {ber Ojcr)^ / \5 l.„-j 

(Pcr[z) ^ - ■ . {z-acr)^e^ , Z -)■ ttcr, 

where arg (;2 — Ocr) E (0,27r), and 

= ^-^^ogz + 0{l), z^O, 


(3.21) 


(3.22) 


where arg^ G (— vr, tt). From the above formula, one can see that lie 0 if t < 0 

and 2 ; approaches the origin such that cos(arg 2 ;) < Figures 3 and In 

particular, we see that is exponentially small as 2 ; —)■ 0, 2 ; G F 2 or 2 ; G F 3 ; cf. 

Figure [^for the contours. In view of (3.19), one can see that the same holds for 


It is worth mentioning that on F 2 and F 3 with \z — 6\ = 6, we have Re “ = 


4 Nonlinear steepest descent analysis 

In this section, we apply the nonlinear steepest descent method developed by Deift 
and Zhou et al. [I2l[l3] to the RH problem for Y. The idea is to obtain, via a series 
of invertible transformations E —)■ T —)■ S' —)■ R, the RH problem for R whose jump 
matrices are close to the identity ones. 
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Figure 3: The shaded region is the region where Re < 0. The left and right pictures 
correspond to cases when t < 0 and t > 0 , respectively. 



Figure 4: The shaded region is the region where Kecpcriz) < 0. Note that this hgure is 
not the exact one for (per defined in (3.15). Here we have rescaled the hgure, especially 
near acr, for better illustration: because the exact value of Ocr is too small as compared 
with ber] see (1.33). 


4.1 The first transformation Y ^ T: Normalization at infinity 


We make use of the ^^-function dehned in (3.13) to normalize the RH problem for Y in 
Section 2.1 when k = n. As g{z) ~ logz for large z, we introduce the hrst transformation 
y —)■ T as follows: 


T{z) = (4.1) 

where I is the Lagrange multiplier in ( |3.16[ ). Then, T solves the following RH problem. 
(Tl) T{z) is analytic in C\F; see Figure [^for F = Fi UF 2 UF 3 ; 
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(T2) The jump condition is 

^n(g-(z)-g+(z)) ^ ,^n{-Vt(z)+g+{z)+g-{z)-l) 


r+(z) = T.(z) 


0 


zn{g+{z)-g- (z)) 


(4.2) 


for 2 ) G Tj, j = 1, 2, 3, where Vt(z) is dehned in (1.7), and Ci = 1, C 2 = a, C 3 = 1 —a; 


(T3) The asymptotic behavior of T{z) at inhnity is 
T{z) = I + 0{l/z) as —)■ cxD. 


(4.3) 


Appealing to the properties of g{z) and in (3.16) and (3.18), the jump matrices 


in (4.2) can be expressed in terms of the function (j)t{z) as follows: 

z eTj \ {acr,b), j = 1,2,3; 


TUz) = Tjz) 


1 

0 1 


^2n{<f>t)^iz) 

0 


(4.4) 


^2n{(j}t)_{z) 


, Z G (flcri b '). 


4.2 The second transformation T ^ S: Contour deformation 

Since {(j)t)± (z) are purely imaginary on (ocr, b), the jump matrix for T{z) on z E {ttcr, b) 
possesses highly oscillatory diagonal entries. To remove the oscillation, we open the lens 
near {acr,b) and introduce the second transformation: 


f Tiz), 


for 2 ; outside the lens shaped region; 


= T{z) ^ j , for z in the upper lens region; 

T{z) ^ 2 n(pt{z) 1 ^ ’ for 2 ; in the lower lens region. 

Then S{z) solves the RH problem 
(SI) S{z) is analytic in C \ S 5 ; see Figure [^for the contours; 


(4.5) 
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Figure 5: Contour S 5 = (ocr, b) U (6, 00 ) Tj- The shaded region is the region where 
Re0cr(2:) < 0. 


(S2) The jump conditions are 




1 0 

^2n4>t{z) ^ 

0 1 
-1 0 
1 

0 1 
1 (1 - 
0 1 

^ ^-2n<t>t(z) 

0 1 


(S3) The asymptotic behavior at inhnity is 
S{z) = I + 0(1/z), as z ^ 00 . 


z G r 4 U Ts, 
Z G (dcr; ) 

^ £ r2, 

) ^ £ r3, 

z G {b, + 00 ). 


(4.6) 


(4.7) 


To study the asymptotic behavior of S{z) for large n, we may exam the signs of 
Ke(j)t{z), to see if the jumps are of the form I plus exponentially small terms. Special 
attention should be paid in the present case since we are dealing with the signed measure 


(3.9). Fortunately, when n is large and t — ter is small enough, we can still determine 


the signs of Re(j)t{z) near the endpoint acr- Similar discussions can be found in 
where modihed equilibrium problems are also addressed. 
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Proposition 4. Let U be a neighbourhood of Ocr- Then, for any £ > 0, there exists a 
6t > 0 such that for all t E R. with \t — tcr\ < St, we have Recftiz) < —e on the upper 
and lower lips of the lens on the outside of U, namely, z G {r 4 U Fs} \ t/. Moreover, 
there exists a positive r, such that Recftiz) > e on {r 2 U Fs} \ U and on [b + r, cxd). 


Proof. In view of (3.12), we see that the factor + d^x + in (3.9) possesses a pair of 
zeros 


x± 




1/4 

(t - + 0{t- ter)- 


One can choose 6t small enough, so that x± G U. Similar to those conducted in 
Prop. 4.2] and [HI P-23], we can prove that the jumps on the portions of F 4 and F 5 , 
outside of U and keeping a distance from the soft edge b, are of the form I plus an 
exponentially small term. 

Next, we estimate Reeptiz) on F 2 UF 3 U( 6 , 00 ) in a straightforward manner, using the 


explicit representations of the 0-functions (3.14), (3.15) and (3.19). In view of (3.19), 


we need only check the critical case for (p^r, and it is readily seen from (3.15) that 


Re {(pcr)+ (a;) = - 


1 r {S - acr)^/‘^{s - bcrY^'^ 


ds > 0 for X > br 


For 2 ; on the semi-circle F 2 ; cf. Figure we take the integration path to be the arc from 
Ocr to z, and adapt the parametrization s = 5 -1- (5e*^, where 5 = acr/2 and 6 G [0, tt). As 
a result, we have 


(peri^z'^ 




(^COS0 (4.8) 


where 9^ G [ 0 , 7 r) such that z = 6 + What is more, we have arg(s — ber) G 

[tt — arcsinA, tt] C [ 37 r/ 4 , tt] on F 2 , where A = 


‘2.hcr etc 


so that arcsinA 0 . 002 l 7 r. 
Hence the argument of the integrand lies in the interval [27r — dvr/S, 27r -|- 7 r/ 2 ), so long 

(ocr) = 0 , and Re(pcr{z) is strictly 
r. The same result holds for z G F 3 . 
It is worth mentioning that Re(pt{z) —>■ -fex) a.s z ^ 0, z E F 2 UF 3 ; see the formula (3.22) 
and the discussion that follows. We note that the estimates on the lens boundaries F 4 
and F 5 can also be obtained from the representations (3.14), (3.15) and (3.19). □ 


as 6 * G [ 0 , 7 r). Therefore, from (4.8) we see that <pcr 
monotonically increasing as z G F 2 goes away from Oc 


4.3 The global parametrix 

Having had Proposition]^ we see from ( |4.6 ) that the jump matrix for S is the identity 
matrix, plus an exponentially small term, for hxed 2 ; bounded away from the interval 
{acr,b). Neglecting the exponentially small terms, we arrive at an approximating RH 
problem for N{z) as follows: 
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(Nl) N{z) is analytic in C\[acr,&]; 

(N2) 

N+ix) = N-ix) \ g ) X e (a^ryff); (4.9) 


(N3) 


N(z) = I + 0(l/z) as z—>■ oo. (4-10) 

The solntion to the above RH problem is constrncted explicitly as 

N(z) = M-^p(z)-^^M, (4.11) 

/ \ 1/4 

where M = (I + iai)/\/2 and q{z) = with arg(2; — Ucr) G (— tt, tt) and 

arg(2; — 6 ) G (—vr, tt). 

The jnmp matrices of S{z)N{z)~^ are not nniformly close to the identical matrix I 
near the endpoints acr and b, thns local parametrices have to be constrncted in neigh¬ 
borhoods of these endpoints. 

4.4 The local parametrix at the soft edge 

The local parametrix at the right endpoint z = b is the same as that of the Lagnerre 
polynomials at the soft edge. More precisely, the parametrix is to be constrncted in 
U{b,r) = {z \ \z — b\ < r}, r being a hxed positive nnmber, snch that 

(a) P^^\z) is analytic in U{b,r)\T,s', see Fignrej^for the contonr S 5 ; 

(b) In U{b,r), P^^^z) satishes the same jnmp conditions as S{z) does; cf. 

(c) P^^\z) fulhls the following matching condition on dU{b,r): 

P^^\z)N-\z) = I + O . (4.12) 

The parametrix can be constrncted, ont of the Airy fnnction and its derivative, as in 
[281 (3.74)]; see also [II1II3]. 

4.5 Local parametrix at the critical point z = Qcr and Painleve 
I 

Now we focns on the constrnction of the parametrix at the endpoint acr- We seek a 
parametrix in U{acr,r) = {z \ \z — acr\ < t}, r being a hxed positive nnmber, snch that 
the following RH problem is satished: 


(4.6); 
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(a) 

(b) 

(c) 


pi^)(z) is analytic in U{acr,r)\T,s; cf. Fignre|^ 

In U{acr,r), P^°'\z) satisfies the same jnmp conditions as S{z) does; cf. (4.6); 
p(“)( 2 ;) fulfils the following matching condition on dU{acr,r): 

p(^\z)N-\z) = I+ 0 asn^oo. (4.13) 


First we make a transformation to convert all the jnmps of the RH problem for (z) 
to constant jnmps. Let ns dehne 

p(a)(^) _ ^ g U{acr,r)\Es, (4.14) 

then it is readily verihed that (z) satishes a RH problem as follows: 

(a) P(“)(z) is analytic in U{acr,r)\T,s; 

(b) In U{acr,r), P^°'\z) satishes the jnmp conditions 

^ £ F 2 , 

^ ^ hs, 

(4.15) 

Z G (cicr 7 ^cr 7 

G 1^4 u r 5 . 


P^Z\z) = pZ\z) 



We are now in a position to constrnct a solntion for the above RH problem by nsing 
the \k-fnnction associated with the Painleve I eqnation, as introdnced in Section 11.11 We 
note that P^°‘\—z) shares the same jnmps with the 4/-fnnction. Then, we establish 
the following conformal mapping: 


fiz) = ( -^(Pcriz) 


2/5 


for 2; G U{acr, t), 


(4.16) 


and r being snfhciently small; cf. (3.21). Indeed, in view of (3.15), one can improve 
(3.21) to obtain 


f (^z') (&cr ®cr) ® (2ncr) ^ i.Z ®cr) ( ^ T dc- 


as z —)■ Qcr- Moreover, we have 


17 f 1 2 , ^ ^ 

25l2(a„-M a. < 


( 4 . 17 ) 




(4.18) 
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where the fractional power takes the principle branch. We also define 


q{z) 


(t ter) 


0o(2:) 


(4.19) 


where the square root takes the principal branch again. It follows from the dehnitions 
of (1)q{z) and f{z) in (3.20) and (4.16) that q{z) is analytic in a neighborhood of z = Qcr 
and 


q{Ojcr) ^ (^cr^cr) ^{per ^cr)^it ter)- 

Moreover, we have 

2 4 

0{n5f{z),n5q{z)) = -nptiz), 


(4.20) 


(4.21) 


where 0{C,s) is the function defined in (1.23). 

With all these preparations, the parametrix near the endpoint Ocr can be constructed 
explicitly as 


P^'^\z) = E{z)'^ (n-^f{z),nlq{z)'^ ^ 

where E[z) is dehned as 

E{z) = N{z)e~^ ^ ^ (n 5 /( 2 :))~ 4 '^ 3 . 

V2 


(4.22) 


(4.23) 


It is ready to see that E{z) is analytic in U{acr,r) and the function P^°''>{z) in (4.22) 
indeed satisfies the matching condition ( |4.13[ ). 

Remark 4. As we have discussed in Section |T3 the solution \k((C;s) exists if and only 
if s is not a pole of ya{s)- Then, to make P^^'pz) in (4.22) well-dehned, we choose 
t — ter = as n —)■ cx), and require nsq^z) is not a pole of |/q,(s). Moreover, from 

the large-C behavior of T((C;s) in ( 1 . 22 ), we can verify the desired matching condition 
on dU{acr,T) in (4.13). 


4.6 The final transformation S ^ R 

With all the parametrices constructed, let us introduce the final transformation 

{ S{z)N-\z), zeC\ {U{a, r) U U{b, r) U S 5 } ; 

^(;,)(p{a))-i(;,), zeU{a,r)\Es-, (4.24) 

S{z){P0’p-\z), zeU{b,r)\Es. 

Then R{z) satishes a RH problem as follows: 

(Rl) R{z) is analytic in C\Sij; see Figure]^ 
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(R2) R{z) satisfies the jump conditions 


(4.25) 


(4.26) 



R+{z) = R-{z)Jr{z), z e T,r, 

where 

{ p(<^)(z)N-\z), zedU{acr,r), 

p(f>)(^z)N-\z), zedU{b,r), 

N{z)Js{z)N-\z), Er \ d{U{a,r, r) U U{b, r)); 

(R3) R{z) satishes the following behavior at inhnity: 

R{z) =1 + 0 {1/z) , as 2 ; —)■ cxD. 


Figure 6: Contour Y^r 


Based on the matching conditions (4.12), (4.13) and the properties of the function 
bt{z) stated in Propositionwe have the following estimates: 


I + O {n , zedU{acr,r), 
jRi^)={ ^ + zedU{b,r), 

I + 0(e-‘='^), zeYrX {dU{a, r) U dU{b, r)), 


(4.27) 


where c is a positive constant, and the error term is uniform for z on the corresponding 
contours. Here we also require that t — ter = 0(n“^/®); see Remark]^ Then, applying 
the standard argument of integral operator and using the technique of deformation of 
contours, see for example nn, we can see that R{z) exits when n is large enough and t 
is close to ter- Moreover, we have 


R{z) =/ + 0(n-^/®). 


(4.28) 


uniformly for z in the whole complex plane. 

This completes the nonlinear steepest descent analysis. 





5 Proof of the main results 

To prove our main results, Theorems and we need more refined asymptotic approx¬ 


imations for R{z) than the one we get in (4.28). 


5.1 Asymptotic approximation of R{z 

For this purpose, we derive an asymptotic approximation for R{z). To simplify the 
statement of our result, we denote 


0 1 
0 0 


0 0 
1 0 


= I 1 and (T_ = 

Proposition 5. Let M = + iai) and R = R{nR^q{z)) ; see (|1.28[). For n —)■ oo, 


t —>■ ter 0'S in ( |1.34[ ), and s* is not a pole of the tronquee solution yais), we have 

1 


R{z) =I + M 


V ui/5 


+ 


2/5 


+ 0 


3/5 


n- 


where 


R^^\z) = 


riR{s* 


z - Or, 




M, 


z e C \ U{acr,r), 


(5.1) 


iR {iRg^ rfR^s*) . 

=(T+ -h ^ -a_, ^ e U{acr,r) 




V7 


Z - Or 


(5.2) 


and 


R^‘^\z) = 


r2 


iVais*) - R{s*))a3, z e C\U{acr,r). 


Z — Ocr 

The constants ri and r 2 in the above formulas are given as 

4/5 


(5.3) 


‘^OjQribcr Oer') 


r2 = 


{2{ber-aer)yR^ 

and S* = nR^q{acr) = — t){2acr)~^ (OerberR^ {ber — Ocr)K of (1.34). 

Proof. From ( |1.22| ) and (4.27), we derive an asymptotic expansion for the jump Jr{z 


(5.4) 


J„(z) = P<“)(.-)iV-‘(.-) = / + 


Jr\V 


k=l 


n 


fc/5 ’ 


2 ; G dU{acr, r), 


(5.5) 


where 


4")(;2) = iV(z)e™^/^T_fc(n5g(z))e-™^/"iV(^)-i(/(^))-'^/^ A; = l,2,--- . (5.6) 

For 2 ) G dU{b,r), Jr{z) is expanded in terms of n~^ for non-negative integers k. Thus, 


we have (5.1) for large n. 


To derive the explicit formulas of R^^\z) and R^‘^\z)., we combine (4.25), (5.1) with 
(5.5). Then one can see that R^R{z) and R^‘^\z) satisfy RH problems as follows: 
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(5.7) 


(i) and fiW(a) are analytic for z G C \ dU (ocr, r); 

(ii) R^^\z) and R^‘^\z) satisfy the jnmp conditions 

R+\z) - R^y{z) = ^^0-+ - z e dU{acr,r) 

qWj Vi 

and 




tH 

77f 




iHg^ 

~7T 


■CT_ 


^ G dU{acr, r), 


(5.8) 


(iii) As z ^ oo, both R^^\z) and R'^‘^'>{z) are of order 0{z ^). 


It follows from the Plemelj formnla that and R^‘^\z) can be represented as 

the Canchy-type integrals on the circnlar contour dU{acr,f') of the right-hand terms in 


(5.7) and (5.8), respectively. The above RH problems can then be solved by conducting 


residue calculations of the Cauchy-type integrals. As a direct consequence we obtain 
(|A^ and ([5)3|. □ 


5.2 Proof of the main theorems 

Now we are ready to derive the large-n asymptotics for the recurrence coefficients and 
the Hankel determinant, as stated in our main theorems. Theorems and 

Proof of Theorems and Tracing back the transformations R ^ S ^ T ^ Y in 
( [T5| ) and ( |4.24D , we have 


Y{z) = (5.9) 

for z close to the origin. To apply the differential identities (2.8) in Lemma we need 


to extract the asymptotics of R(0) when fc = n, n —)■ cx), t —?• ter as in (1.34), and s* is 


not a pole of the tronquee solution ?/a(s). From the explicit formula of N{z) in (4.11) 


and the approximation of R{z) in (5.1), we have 


2 \ i[d- 2 ) d+ 2 ^ 


riR{s* 


d' d 

/2 / 


r 2 iVais*) - n\s*)) 

2a.r 77,2/5 


2acrdn^/^ 

\-d ^{d+\) 

\-d 


—i 1 
1 i 


+ 0 


3/5 


w 


(5.10) 


/ X 1/4 

where d = and the constants rj are dehned in (5.4). Now (2.8) in Lemma 

implies that 


= -A{R{0)N{0))u{R{0)N{0)hi. 


(5.11) 
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Then the asymptotic formula (1.35) for the Hankel determinant follows immediately 


from the identities (5.10) and (5.11). This completes the proof of Theorem]^ 

To derive the asymptotics of the recurrence coefficient and the leading coefficient 
7n,n, we use the relations 




21 , 


7^ = 

I n.n 


2Tii(Y_ 


(5.12) 


1)12 


where Tli is the coefficient of the 0{1/z) term in the asymptotic expansion of Y{z)z 
that is, 


Y{z)z 


—nos _ 


OO 


i=i 


z^ 


as z ^ OO. 


(5.13) 


Therefore, we also need to expand N{z) and R{z) as z ^ oo. Again, using the explicit 


formula of N{z) in (4.11) and the asymptotic approximation of R{z) in (5.1), we have 

N(z) = / H- - + O ( ^ , with N_i = zib — acr)M~^(T^M, (5-14) 

z \z J 4 

and 


R{z) _ / + ^ + 0 ) , 


where 


-1 ^1/5 + ^2/5 + WTsi 


(5.15) 


(5.16) 


Recalling the values of ri and r 2 in (5.4), and the identities M = 4 


0 i 


and M = ( ), we have from (|5.9|) and (|5.12|) that 

—i 0 


—i 1 
1 i 


f^n,n — {R-1 + A^-i)i 2(-R-1 + A^-l)21 

{bcr ®cr) (2.(lcr(bcr ®cr)) ^ Vai^S ) 1 


16 


n 


2/5 


+ 0 


^ 3/5 


This is ( |1.37[ ). 

Similarly, the asymptotics for the leading coefficient can be obtained. Indeed, we 
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have 


^—nl 


^2 = 
/ n,n 


27ri {R-i + A^-i)i2 

^brr-arr. r{H{s*) ^r 2 (|/a(s*) - H(s*)^) 


2ni 

2e-”' 


'cr . , 

-% + 


4 " ' 2nh5 

2ir{K{s* 


+ 


'^{bcT ^cr) 
26-"-^ 


+ 


^2/5 

4r2(2/a(5*)-?{(s*)2) 


+ 0 


^3/5 


-1 


1 + 


(6cr - acr)nR^ {b, 

2ir{H{s*) 


cr ^cr 


)n2/5 


+ 0 


n3/5 


'^(bcr V {per 

4r^-H(g*)^ + 4r,{ye,is*) - U{s*f){ber - ae 
Per - aerYn^/^ 


+ 0 


n 


3/5 


(5.17) 


which gives us (1.38). 


Finally, we derive the asymptotics for an,n- Note that, from (2.6), we have 
an.n = 2ir!i7/„Fn(0)Fi2(0). 


(5.18) 


To hnd the asymptotic behavior of Fii(0)Fi2(0), we have from (5.9) and (5.10) that 

ni(o)yi2(o) 

= e’^'(i?(0)iV(0))n(7?(0)iV(0))i2 

nl 


= e 


—i 


= —e 


--((Y - rf-2)_ nUjs*) 

rln{s*Y ^ r 2 {ya{s*) -'H{s*Y)per + acr) '\ 1 ^ ^ / 1 


AOjerY ^crbe7 
nl per ®cr)^ 

4'\/ (Xerber 

rl'H{s*Y 

^erPer ^cr') 


1 - 

+ 


20‘erY Berber 
2iri'H{s*) 
per ^cr)^ 

‘2r2{yap*) - Hp*Y)Per + 


^2/5 


^3/5 


^erPer ^cr) 


^2/5 


+ 0 


n3/5 ) I ■ 


Then, the asymptotic approximation for an,n in (1.36) follows from the above formula 

□ 


and (15.171). 


This completes the proof of Theorem 
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